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Abstract 

A re-examination of the atmospheric cascade profile of the highest energy cosmic 
ray is presented. The study includes air-shower simulations considering different 
cross sections, particle multiplicity and variation of the hadronic-event-generator to 
model interactions above 200 GeV. The analysis provides evidence that a medium 
mass nucleus primary reproduces the shower profile quite well. This result does not 
support the idea, increasingly popular at present, that the highest energy particles 
are protons, derived from the decay of supermassive relic particles. On the other 
hand, we show that debris of relativistic super-heavy nuclei, which can survive a 
100 Mpc journey through the primeval radiation are likely to generate such a kind 
of cascade. 
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Shortly after the microwave background radiation (MBR) was discovered [1], it 
became clear that ultra high energy cosmic rays (CRs) undergo reactions with 
the relic photons yielding a steep drop in their energy attenuation length. 
Therefore, if the CR sources are all at cosmological distances, the energy 
spectrum would exhibit a Greisen-Zatsepin-Kuz'min (GZK) cut off slightly 
above 100 EeV (1 EeV = 10 18 eV) [2]. The spectral cut off is less sharp for 
nearby sources, viz., distances < 100 Mpc. 

Over the last few years, several giant air showers have been detected, which 
imply the arrival of "super-GZK" particles [3]. Namely, the Haverah Park 
experiment reported several events with energies near or slightly above 100 
EeV [4]. On October 15, 1991 the Fly's Eye experiment recorded the highest 
energy cosmic ray ever detected on Earth (E ~ 300 EeV) [5] Q More recently, 
the Akeno Giant Air Shower Array (AGASA) has reported several remarkable 
CR events scattered across half the sky dramatically confirming that the CR 
spectrum does not end with the GZK cut off [7]. 

Many theories have been suggested to explain the origin and nature of these 
particles. At the most general level, the production mechanisms can be classi- 
fied into two distinct groups: "bottom-up" acceleration arena, and "top-down" 
decay scenario, each of them branching out into various classes. 

The "bottom-up" mechanism usually implies the stretching of rather well 
known acceleration processes to their theoretical limits (and sometimes be- 
yond). It involves particle acceleration in the accretion flows of cosmological 
structures. Examples are one-shot acceleration in active galactic nuclei (AGNs) 
[8], large scale shocks resulting from structure formation [9], acceleration in 
quasar remnant cores of nearby giant elliptical galaxies [10], collisions of galax- 
ies [11], acceleration in catastrophic events associated with gamma ray bursts 
[12], and acceleration in hot spots of extragalactic radio sources [13]. 

On the other hand, the "top-down" scenario (mischievously) escapes from the 
acceleration problem by assuming that charged and neutral primaries simply 
arise by quantum mechanical decay of supermassive elementary X particles 
{mx ~ 10 25 — 10 28 eV) [14]. Sources of these exotic particles at present could 
be topological defects (TDs) left over from early universe phase transitions 
associated with spontaneous symmetry breaking underlying unified models of 
high energy interactions. Due to their topological stability, the TDs (magnetic 
monopoles, cosmic strings, domain walls, etc.) can survive forever with X 
particles trapped inside them. Nevertheless, from time to time, some TDs can 
be destroyed through collapse, annihilation, or other processes, and the energy 
stored is released in the form of massive quanta that would typically decay 
into quarks and leptons. The subsequent hadronization of the quarks would 

1 Near the arrival direction of this event the Yakutsk experiment detected an ex- 
tensive air shower whose energy may also reach 300 EeV [6]. 
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produce jets of hadrons containing mainly light mesons (pions) together with a 
small fraction (3%) of nucleons. Thus, the "top-down" injection spectra would 
be dominated by gamma rays and neutrinos originated by pion decays. 

There have been arguments for some time that the cosmic photon background 
makes the universe opaque to ultra high energy gamma-rays [15]. The most 
efficient target for these rays to produce e + e~ pairs are background photons 
of energy ~ m 2 e / E < 10 _6 eV ~ 100 MHz. The mean free path for a 100 EeV 
(1000 EeV) photon to annihilate on the radio background is believed to be 
10-40 Mpc (1-5 Mpc) [16,17]. Consequently, the measurement of the diffuse 
gamma-ray background imposes several constraints on the total amount of 
electromagnetic energy density injected and recycled by cascading to lower 
energies, yielding bounds on TD injection spectra [18]. It is also notewor- 
thy that the expected high diffuse neutrino flux above 1 EeV sets additional 
bounds on the "top-down" decay scenario [19]. 

The identity of the primary particle certainly deepens the mystery. Strictly 
speaking, a 100 EeV proton traverses a nearly straight line through the galactic 
or extragalactic magnetic fields and yet no compelling local sources have been 
identified within the expected scattering cone of the highest energy event [20]. 
Furthermore, the atmospheric shower profile is reported by Fly's Eye [5] to 
roughly fit that of a proton, or a gamma-ray that peaks somewhat late in the 
atmosphere [16]. It is also reported that a medium mass nucleus seems to fit 
this shower profile. 

Heavy nuclei have their own merits. They can be accelerated to higher terminal 
energies because of their larger charge. Additionally, a heavy nucleus would 
have a larger scattering cone since they can be deflected more by the magnetic 
fields on the way to the observer. However, there are still other problems 
left in interpreting all the "super-GZK" events as heavy nuclei: It has long 
been known that above 200 EeV nuclei should be photodissociated by the 
MBR within a few Mpc [21], and there are few nucleus-emitting-sources in 
the neighborhood of the Earth. More speculative primaries have also been 
proposed, such as dust grains [22] or even neutrinos [16]. On the one hand, the 
grains turn out to be unstable with respect to the development of a fracture. 
The path length up to the first break-up turns out to be much less than the 
characteristic size of the Milky Way [23]. In addition, the depth of shower 
maxima registered by the Haverah Park and Volcano Ranch experiments do 
not correspond to the expected picture of showers initiated by relativistic 
specks of dust [24]. Neutrinos, on the other hand, can traverse unscathed 
through the primeval radiation, but they would predominantly give rise to 
deeply penetrating showers in the atmosphere, again, in disagreement with 
observations. Putting all these together one can naively say that aside from 
some nucleus species, the broad phylum of particles does not seem to fit the 
cascade development of the highest energy event. 
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All in all, the determination of the chemical composition appears crucial in 
discriminating between the possible origin of the highest energy cosmic rays. 
In light of this, we decided to re-examine the shower profile reported by the 
Fly's Eye experiment. In this letter we present an extension of the previous 
analysis [16] studying the effects of different atmospheric-cascade-models and 
parameter values. 

The Fly's Eye observes an air shower as a nitrogen fluorescence light source 
which moves at the speed of light along the path of a high energy particle 
traversing the atmosphere [25]. In other words, it directly detects the longitu- 
dinal development of the air cascade. The simulation of the shower evolution 
depends sensitively on the first few interactions, necessarily related to the 
quality of our "understanding" of hadronic collisions. Soft multiparticle pro- 
duction with small transverse momenta with respect to the collision axis is 
a dominant feature of most events in high energy hadronic interaction. De- 
spite the fact that strict calculations based on ordinary QCD perturbation 
theory are not feasible, some phenomenological approaches successfully take 
into account the main properties of soft diffractive processes [26-28]. In these 
models, the interactions are no longer described by single particle exchange, 
but by highly complicated collective modes known as Reggeons [29]. The slow 
growth of the cross section with the centre of mass energy requires a dom- 
inant contribution of a special Reggeon, the Pomeron [30]. As a result of 
the diffractive interactions, constitutents get excited and produce particles, 
which, in the fragmentation phase, are modelled by strings [31]. The physical 
picture at super-high energies is quite different. Interactions are not fully de- 
scribed within the framework of standard Gribov-Regge theory [29]. At these 
extremely high energies, apart from the usual excitations of the participants, 
hard parton-parton collisions (between the beam and target hadrons) take 
place. A suitable scenario for the joint description of both, soft and semi-hard 
hadronic physics has yet to see the light of day. The well known codes QGSJET 
[32], sibyll [33], and dpmjet [34] represent three of the best simulation tools 
to model hadronic interactions at the highest energies. Let us now concentrate 
on the first two programs, and a discussion on DPMJET will follow.0 

The underlying idea behind sibyll is that the increase in the cross section 
is driven by the production of minijets [33]. The probability distribution for 
obtaining N jet pairs (with > p™ m , being p™ in a sharp threshold on 
the transverse momentum below which hard interactions are neglected) in a 
collision at energy ^/i is computed regarding elastic pp or pp scattering as 
a difractive shadow scattering associated with inelastic processes. The algo- 
rithms are tuned to reproduce the central and fragmentation regions data up 



2 Improvements in the hadronic models are also under way. NEXUS, a joint enterprise 
by the authors of venus [28] and QGSJET is expected to become available soon. The 
preliminary nexus skeleton has already been described in [35]. 
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to pp collider energies, and with no further adjustments they are extrapo- 
lated several orders of magnitude. On the other hand, in qgsjet the theory 
is formulated entirely in terms of Pomeron exchanges. The basic idea is to 
replace the soft Pomeron by a so-called "semihard Pomeron" , which is defined 
to be an ordinary soft Pomeron with the middle piece replaced by a QCD par- 
ton ladder. Thus, minijets will emerge as a part of the "semihard Pomeron" , 
which is itself the controlling mechanism for the whole interaction [32]. The 
different approaches used in both codes to model the underlying physics show 
clear differences in multiplicity predictions which increase with rising energy 
[36]. It is important to stress that although the differences seem to become 
washed out as the shower front gets closer to the ground, the footprints of 
the first hadronic collisions are still present in the longitudinal development. 
Consequently, the first generation of particles plays a paramount role in deter- 
mining the chemical composition of the "super-GZK" event reported by Fly's 
Eye. 

In order to determine the progenitor of the highest energy event, we carried out 
a Monte Carlo simulation of the atmospheric shower profile using the aires 
program (version 2.2.1) [37]. The simulation includes interaction with both the 
Earth magnetic field and nuclei in the atmosphere. We have taken into account 
possible representative variations in the transport code, cross sections, particle 
multiplicity and hadronic models by using SIBYLL and QGSJET to generate 
interactions above 200 GeV. 

We generated several sets of nucleus showers with the mass equally spread in 
the aires range (up to 56 Fe). In addition, air-showers of superheavy nuclei 
were simulated following the superposition model. This model assumes that 
an average shower produced by a nucleus with energy E and mass number A is 
indistinguishable from a superposition of A proton showers, each with energy 
E/A. The particles were injected at the top of the atmosphere (100 km.a.s.l) 
with zenith angle 43.9°, and the magnetic field was set to reproduce that 
prevailing upon Fly's Eye experiment [38]. All shower particles with energies 
above the following thresholds were tracked: 750 keV for gammas, 900 keV for 
electrons and positrons, 10 MeV for muons, 60 MeV for mesons and 120 MeV 
for nucleons and nuclei. The results of these simulations were processed with 
the help of the aires analysis programs. 

The data analysis is performed by means of a y 2 test [39]. We assume that 
the set of measured values by Fly's Eye are uncorrected, i.e, any depth mea- 
surement is independent of each other. Then, we make use of the quantity 

N I |2 

where N is the total number of points in the analysis, a Xj is the error on 
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the Xjth coordinate, Xj is the measured value of the coordinate, and aj the 
(hypothetical) true value of the coordinate. Generally speaking, the simulated 
shower profile would give a good representation of the data if the \ 2 is roughly 
equal to the number of data points, i.e., x 2 ~ 12. For a gamma-ray primary 
we obtain Xsibyll = 1235.82, and Xqgsjet = 1175.39. The results for different 
primary nuclei are shown in Fig. 1. 

Definite conclusions on the chemical composition cannot be reached, mainly 
because of large fluctuations from model to model. Without belaboring the 
point, the main differences between the hadronic-event-generators can be as- 
cribed to two different parameters: the inelastic cross section, and the energy 
fraction of the parent particle converted into the production of secondary par- 
ticles, the so-called "inelasticity". In order to match collider data, the brutal 
truncation of the soft hadronic processes (sibyll) enforces interactions with 
a large average fraction of the energy going into the leading baryon (elastic- 
ity) and a small fraction into secondaries (inelasticity) than those governed by 
QGSJET. As a consequence, sibyll's showers require more generations of par- 
ticles undergoing hadronic collisions, yielding a delay in the electromagnetic 
shower development (strongly correlated with 7r° decays). However, both the 
cross section, and the inelasticity seem to affect the shower development in 
a similar fashion since increasing inelasticity as well as cross section lead to 
faster cascades. Therefore, primaries with very different combination of cross 
section and inelasticity are able to describe the given longitudinal develop- 
ment. In particular, the large elasticity of sibyll yields a minimum on the 
X 2 above the aires nucleus range. Here, the reader should keep in mind the 
crudeness of the superposition model for air shower simulations. In the simu- 
lations with QGSJET+ aires, the x 2 reaches a minimum at 41 Ca, leading to a 
plain representation of the shower profile (see Fig. 2). 

The total error in the energy determination (systematic and statistical uncer- 
tainties added in quadrature) reported by Fly's Eye is 93 EeV [5]. The next 
step in the analysis is to check whether the above results could be improved 
for any particular energy (within the error) by forcing the primary to be a 
proton. To this end we generate several sets of showers, varying the primary 
energy within the error box of the event. The x 2 always remains very large. 
For example, in Fig. 1 we show the result (as a star) obtained for the upper 
limit of the energy for comparison. The last exercise is to analyze the sensi- 
tivity of the shower simulations. To handle shower-to-shower fluctuations we 
separately analyze each shower initiated by a proton: after 50 showers all the 
X 2 s are above 20. Putting all this together, the point to be made is this: Even 
though the primary chemical composition remains hidden by the hadronic 
interaction model, it is evident that in both cases the shower is completely 
inconsistent with a proton, or a gamma-ray primary. 

At this stage, it is important to stress that although the latest version of dpm- 
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JET was improved for operation up to primary energies of 10 21 eV (per nucleon 
in the lab. frame) [34], it has not yet been effectively implemented into the 
standard transport codes [40]. However, one does not expect strong deviations 
in the longitudinal development if the hadronic collisions are modelled with 
dpmjet [41,42]; the results for the x 2 should lie somewhere between those of 
aires+SIBYLL and aires+QGSJET (see, in particular Fig. 5 of Ref. [41]). 



Let us end by reviewing the current observational status, and discussing the 
nature of the highest energy CRs. 



• Up to now, the AGASA experiment reported 7 events above 100 EeV. Large 
shower to shower fluctuations make it difficult to get a definite dependence 
of composition from muon density at 1000m from the shower core (see Fig. 
3 of Ref. [43]). The situation is further complicated if we take into account 
that not only two mass components (iron-proton) are allowed, but that 
many A- values could be present. There is no candidate shower of gamma- 
ray primary in this sample [43]. 

• The density profile of the highest energy event detected by the Yakutsk ex- 
periment is almost dominated by muons, which argues against proton/gamma- 
ray initiation [6]. 

• The atmospheric profile of the highest energy cosmic ray is inconsistent 
with that initiated by a proton or a gamma-ray. Moreover, an analysis of 
the frequency distribution of the depth maximum values recorded by Fly's 
Eye, indicates that in the energy range 3-10 EeV, there is a significant 
fraction of nuclei with charge greater than one [44]. 

• Very recently, the Haverah Park experiment reported data from inclined 
showers [45]. At the "super-GZK" energies, the resulting spectrum is con- 
sistent with nucleus primaries (see Fig. 2 of Ref. [45]). We stress once more 
that above the photodisintegration threshold different type of debris may 
reach the atmosphere; consequently, the primary mixture is far from just 
protons and iron nuclei. This sample also states evidence against the photon 
hypothesis. 



Whatever the source (s) of the highest energy cosmic rays, because of their 
interactions with magnetic fields, and also with the universal diffuse radia- 
tion backgrounds permeating the universe, the energy spectra, composition 
and arrival directions are affected by propagation. In particular, in the case of 
heavy nuclei, the composition of the arriving particles may differ substantially 
from that of the source. For instance, an iron nucleus which is injected into 
the intergalactic medium with a Lorentz factor of 7 = 4 x 10 9 , after prop- 
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agating m 100 Mpc will end up with A « 30.0 Such metamorphosis is the 
result of two basic interactions: photodisintegration and pair creation. The 
latter reduces the Lorentz factor significantly during propagation avoiding the 
complete disintegration of the nucleus. Stecker and Salamon have pointed out 
that a heavy nucleus emitted by the strong radio galaxy NGC315, which is 60 
Mpc away and within the scattering cone of the highest energy event detected 
by AGASA, would have an energy cut off ~ 130 EeV, which may be within 
the uncertainty in the energy determination for this event [46]. Furthermore, 
it was already noted that a positively charged particle injected by the next- 
door galaxy M82 may have been deflected towards the arrival direction of 
the highest energy event [48]. Note that although M82 is not a strong radio 
galaxy, it has been described as the archetypal starburst galaxy (supernova 
rate as high as 0.2 - 0.3 yr" 1 ) and as a prototype of superwind galaxies. In 
this picture, the central region of the galaxy harbors a population of heavy 
ions (mainly produced in supernova explotions) that could be re-accelerated 
to ultra high energies at the terminal shocks of galactic superwinds generated 
by the starburst. It is noteworthy that M82 could also be responsible for the 
"super-GZK" Yakutsk event. In addition, it is interesting to point out that 
the regular component of the galactic magnetic field acts as a giant lens for 
nuclei where E/Z < 50 EeV, yielding significant deflections of CR trajectories 
[49]. 

Finally, we would like to stress another remarkable feature of nucleus pri- 
maries. It was already noted that nuclei heavier than iron could be accel- 
erated to super-high energies [50]. Moreover, the dominant mechanism for 
energy losses in the bath of the universal cosmic radiation is the photodisin- 
tegration process, which was studied in detail for the case of 197 Au (see Fig. 
2 of Ref. [50]). In particular, the photodisintegration rate R of a gold nucleus 
with Lorentz factor 7o = 1 x 10 9 ( 7o = 2 x 10 9 ) is R = 3.22 x 10~ 15 s" 1 
(R = 1.55 x 10~ 14 s _1 ), and the energy loss time is defined by 

-i = ±dE = ldr L R 
T E dt ydt A' U 



The energy loss due to photopair production could be estimated to be roughly 
Z 2 /A times higher than that of a proton with the same Lorentz factor, i.e., 
for 70 « 10 9 , T~+ e - = l/jdj/dt ~ 2 x 10~ 18 s _1 [51]. For a gold nucleus then 
we get, rs;6x 10 16 s (r ~ 1 x 10 16 s). Now, using the previous formulae for 
r, it is straightforward to show that a gold nucleus emitted with a Lorentz 
factor « 10 9 , after propagating 100 Mpc still has an energy above 100 EeV. 

3 It should be noted, that the photodisintegration history of the nucleus was per- 
formed using the single/double- nucleon emission energy threshold of ref. [46]. The 
value of A is thus slightly different to that of Ref. [47], because the increase of the 
threshold energy lengthens the propagation distance. 
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Specifically, gold nuclei injected with 197 EeV, will mutate to A pa 166 — 
168 with an energy pa 167 EeV, whereas for E = 394 EeV, the surviving 
fragment would have A = 89 — 91 with an energy pa 180 EeV. Note that both 
values are within 1 standard deviation of the highest energy event recorded 
by AGASA (E = 210 EeV, whose accuracy on energy determination is 30% 
[7]). Furthermore, in this case the deflections of the CR trajectories would 
be enough to conceal the source location. The above valuation together with 
our previous analysis suggest that elements heavier than iron are likely to 
generate some of the "super-GZK" events. It perhaps should be interesting to 
try to simulate showers initiated by super-heavy nuclei considering the internal 
nuclear structure, specially when the hadronic interactions are modelled with 

SIBYLL. 

In summary, nucleus induced showers are in complete agreement with the 
experimental data reported by the groups of Fly's Eye, Yakutsk, AGASA 
and Haverah Park. Moreover, if the highest energy CRs are heavy nuclei, one 
cannot yet rule out that extragalactic/galactic magnetic fields could tangle 
up the particle paths, camouflaging the exact location of sources. The astro- 
physics parameters governing CR propagation are still somewhat uncertain, 
and better statistics on the arrival directions, energy spectra and composition 
are indispensable for understanding the origin of the highest energy particles 
ocurring in nature. Nonetheless, it does not seem reasonable to invoke "top- 
down" models to explain the existing data since they would fail to explain the 
nucleus-component, certainly present at the end of the spectrum. 
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Fig. 1. Behavior of x 2 as a function of number of nucleons A. The white star stands 
for the x 2 of a 355 EeV proton induced shower (aires+QGSJET, see main text). 
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Fig. 2. Atmospheric cascade development of a 300 EeV nucleus induced shower, 
superimposed over the Fly's Eye data. 
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